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ABSTRACT 


Independent Modal Space Control (IMSC) is a technique that is often used for 
the control of large order structural systems. The pertinent optimization 
problem in the simultaneous design and control of structures is a min - min 
problem that minimizes with respect to the structural design variables, the 
minimum value of the performance index with respect to the control forces ob- 
tained using the IMSC technique. The minimization process requires derivatives 
of eigenvalues and eigenvectors with respect to the design variables. These 
derivatives can be computed by a rather involved analytical procedure or a 
relatively simple finite difference procedure. This paper examines the computer 
cost effectiveness of these two procedures for the derivative calculations. 


INTRODUCTION 


One of the objectives of structural control is to suppress undesirable 
motion resulting from some unavoidable excitation such as onboard machinery or 
docking maneuvers. In active control the motion of structure is sensed and 
suitable forces are applied to reduce and ultimately eliminate the undesirable 
motion. In optimal control the forces are applied such that a preselected 
performance index is minimized. The solution of the optimal control problem 
requires the solution of the matrix Ricatti equation. Because of the difficul- 
ties encountered in numerical computations, the solution of the matrix Ricatti 
equation is not feasible for large order systems. For large order systems, an 
alternate method known as the Independent Modal Space Control (IMSC) [1] is more 
suitable. 

In the IMSC method, the control forces are specified in the modal space 
instead of in the physical space. Also by suitably choosing the modal control 
forces, each mode of vibration is controlled independently of the other modes. 
The performance index is assumed to be of the form 


i 



where a is the number of modes controlled and J r is the performance index 
associated with the r-th mode and has the definition 
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u> is the frequency of the rth mode and R >0 is the penalty parameter imposed on 
tne control effort. A higher value of R will result in a smaller control force 
in the modal space and vice versa. r 

The modal coordinates £ , r = l,...t are related to the displacement vector u, 
by the relation 


u n X£ (4) 

where X is the modal matrix, having as its columns the eigenvectors, obtained by 
the solution of the eigenvalue problem 

K X = u) 2 MX (5) 


£ r and n r satisfy the constraint equations 
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£ (t) + w £ (t) = z (t) (6. a) 
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t r (0) = n r (0) u> r 


(6.c) 


where Z Q = X T F (7) 

is the modal control vector. Minimization of J in Eq. (2) with the differential 
constraint equations given by Eqs. (6) leads to a 2 x 2 matrix Ricatti equation 
that can be solved analytically for t^ = 

For this case, the control force is given by [2] 



( 8 ) 
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and the solution of the closed loop modal equations (6) for the controlled modes 
gives 


5„(t) = e (a. cos et + p. sin 0t) 


(9) 


with 


a l ■ 

u) 

r r 

(10) 

0 = 

(•<), 

(11) 

°1 B 

^ro 

(12) 

Pi = 

(“d) (^ ro * r + nr ^0 

(13) 

X r - 

f 22 

2u V 

(14) 

(“<O r = 

f 2 - f„.l/2 

0°r - “r f 21 - -^f) 

(15) 

f 21 = 

R r : k 21 , f 22 = “ R r lk 22 

(16) 

k 12 

■ a) R + ( to 2 R 2 + R 
r r l r r r J 

(17) 


k 22 = ( R r- 2 “r R ? + 2 “r R r (“r R r + R rf* j'* 


(18) 


Substitution of Eqs. (9) through (18) into Eq. (2) followed by its integration 
with t^ = oo yields 


0 ) 
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with 


2 k 12 k 22 , 

2 D ’ J 22 
w r R r 



E = J n a 2 + J 22 a 2 + J 12 aj « 2 

F = J 1]L 0 * + J 22 P 2 + J 12 P 1 P 2 

G “ 2 J U a l P 1 + 2 J 22 a 2 p 2 + °12 ( a l p 2 + a 2 P l) 


( 24 ) 
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Simultaneous Design and Control 


The process of simultaneous design and control of structures is a min - min 
problem that minimizes with respect to the structural design variables the minimum 
value of the performance index with respect to the control forces. 

The minimum value of the performance index with respect to control forces is 
given by Eq. (1) wherein J can be evaluated from Eqs. (8) through (29). In the 
process of minimization of this minimum value of the performance index, its 
derivatives with respect to the design variables are required. These derivatives 
can be evaluated explicitly by a laborious, even though straightforward 
differentiation of Eqs. (8) through (29) with respect to the design variables 
provided the derivatives of the eigenvalues u> and the eigenvectors x , r = 1, 
2...Z are available. The other alternative isto calculate the derivatives of J 
with respect to the design variables by using say the forward difference scheme. 
The latter is easily programmable since no explicit derivatives of the eigenvalues 
and eigenvectors with respect to the design variables are then required. The 
thrust of this paper is a comparison of the computational cost and the efficiency 
of the two procedures for calculating the derivatives of the performance index. 

Before we elaborate on this comparison however, we will digress and discuss 
the calculation of the derivatives of the eigenvalues and eigenvectors using the 
well known Nelson's method [3]. 

Derivatives of Eigenvalues and Eigenvectors with Respect to the Design Variables 

Purely from a computer programming point of view the simplest and the most 
straightforward though not necessarily the most efficient way to compute the 
derivatives of eigenvalues and eigenvectors is by using finite differences in 
particular the forward difference scheme with an appropriate step size [4]. The 
main disadvantage of the forward difference scheme is that it requires the 
solution of an eigenvalue problem once for each design variable. This could be a 
computationally expensive process. Furthermore, to obtain an accurate value of 
the computed derivatives, the eigenvalue problems need to be solved with a high 
degree of precision. 

2 

The eigenvalue u> , and the eigenvector x , of the previous section are obtained by 
the solution of the eigenvalue problem. 

<*y M x r = K x r , r = 1, 2...£ (30) 


where M and K are the assembled mass and stiffness matrices respectively of the 
finite -element model of the structure. The mode shape x is normalized with 
respect to the mass matrix M as r 


x M x = 6 
r s rs 


(31) 


752 



wherein 6 is the Kronecker delta. Differentiating Eqs. £30) and (31) with 
respect to rs a design variable p. for a particular eigenpair (o/,x ) with distinct 
eigenvalues one obtains J 



where use has been made of the symmetry of the mass matrix M. 

To obtain the derivatives of eigenvalues, Eq. (32) is premultiplied by x^ 
followed by the use of Eq. (31) to yield 



dK 

dp. 


U) 


dMA 

d V 


(34) 


To obtain the derivatives of the eigenvector x Eqs. 
bined as ’ 


(32) and (33) are com- 



- \ M 


M 



Equations (35) could be solved for both the eigenvalue and the eigenvector 
derivatives except that the principal minor K - ur is singular. To circumvent 
this apparent difficulty. Nelson [5] proposed a method that temporarily imposes 
the normalization equation (31) by the requirement that the largest component of 
the eigenvector be equal to one. If the re-normalized eigenvector is denoted by 
x and it is assumed that its largest component is the m-th one, then Eq. (31) is 
replaced by 


x 


rm 


= 1 


(36) 
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and Eq. (33) is replaced by 



For x Eq. (32) reduces to 


( 37 ) 
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(38) 


Equation (37) is now used to reduce the order of Eqs. (38) by deleting the m-th 
row and m-th column. When the eigenvalue « is distinct, the reduced system is 
not singular and can be solved by a standarS technique for the derivative _vector 
dx r /dp.. The required vectors x and dx /dp. are then obtained from x and 
dxj,/dp^ by the following easily verified refaticms 
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In finite element computer codes that exploit the sparsity structure of the K 
and M matrices, it may be inconvenient to obtain the re-normalized vector, x by 
setting the largest component to unity. Such a scheme necessitates the recal dila- 
tion of the sparsity structure. Instead, it is more convenient to obtain the 
re-normalized eigenvector, x by setting 


x 


rn 


= 1 


(41) 


where n is the order of the matrices K and M. 

As mentioned previously, the derivative of the eigenvector x with respect to 
the design variable p. can also be calculated by the forward difference scheme 

\J 


dx r ^ r V h ' ^ P J 

35J = h 


(42) 
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where (x 1 is the eigenvector calculated at p... . In order to assess the 

r Vh j+h 

accuracy of the forward difference scheme relative to Nelson's analytical method, 
an error measure is defined as 


* <"> ■ 2 2 2 
j r 1 


where 


f 8x rn 


rfi 
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and 
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•1 - (fill 2 

(43) 
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i the eigenvector derivatives by 

the analyti- 


V^Pj 

cal and the forward difference scheme respectively. The 


error e is summed over 
all the components of the eigenvector, over the mode shapes controlled and over 
all the design variables. 


Application to a Stiffened Composite Plate 

A laminated composite square plate reinforced by two stiffeners placed 
symmetrically with respect to the laminate midplane along the two centerlines of 
the plate is considered. As in reference [6] the plate is discretized using a 
mesh of 8 noded isoparametric, shear deformable plate bending elements. Assuming 
the plate is simply-supported along all its four edges, the resulting finite- 
element model has a respectable (from control engineer's point of view) 127 
degrees of freedom and thirteen design variables consisting of five discrete fiber 
orientations and eight continuous stiffener cross-sectional areas. 

Table 1 provides an assessment of the error e as a function of the step size h 
for the finite difference derivative calculations. As expected, the error 
decreases with a decrease in h and then begins to increase as a result of machine 
roundoff. 

A comparison was made of the computational cost for the calculation of the 
eigenvector derivatives using .Nelson's method and the finite difference scheme. 
Using Nelson's method to compute the gradient of the three eigenvectors with 
respect to the thirteen design variables the required CPU time was 17.2 seconds. 
To compute the eigenvector gradients using forward differences several eigenvalue 
problems need to be solved. Using subspace Iteration in conjunction with the 
Jacobi method [7] for the solution of the eigenvalue problem, the total time for 
the required gradient calculations was 39.5 seconds. Note that the design vector 
has thirteen variables, and it was necessary to solve the perturbed eigenvalue 
problem thirteen times. Since the solution of the unperturbed eigenvalue problem 
provides an excellent guess for the eigenvalue of the perturbed system, an inverse 
iteration scheme [7] in conjunction with shifting of the stiffness matrix K can be 
used to accelerate the solution process. Using such a strategy, the CPU time 
required for the calculation of the eigenvector gradients using forward 
differences was down to 27.6 seconds. 
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Thus, in relation to the analytical method the computational cost of the finite 
difference calculation of the eigenvector gradients is not at all prohibitive. On 
the other hand, in spite of this modestly higher computational cost, the 
simplicity of the calculation of the eigenvector gradients using forward 
difference scheme is overwhelming. However some caution must be exercised when 
using inverse iteration in conjunction with shifting of the stiffness matrix. It 
should be noted that the normalization scheme x 7 M x = 1 fixes only the magnitude 
of the eigenvector and if x is an eigenvector, r then-x is also an eigenvector of 
the system. Hence, when eigenvectors of the perturbed system are computed, care 
must be taken to choose the eigenvector (x ) such that 

Pj +h 


<x p ) 


T 

Y h 



> 0 


This can be done very easily in practice by simply calculating the above dot 
product and changing the sign of the vector x r if the dot product is negative. 

Table 1. Error e as a Function of the Step Size h 


h 

e 

0.156 x 10 -1 

0.79 x 10" 6 

0.781 x 10~ 2 

0.21 x 10“ 6 

0.390 x 10" 2 

0.54 x 10‘ 7 

0 195 x 10“ 2 

0.13 x 10" 7 

0.976 x 10" 3 

0.34 x 10" 8 

0.488 x 10~ 3 

0.98 x 10 _9 

0.244 x 10" 3 

0.80 x 10 -9 

0.122 x 10' 3 

0.25 x 10“ 8 

0.610 x 10‘ 4 

0.101 x 10' 7 

0.305 x 10" 4 

0.404 x 10‘ 7 

0.152 x 10~ 4 

0.162 x 10 _6 
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Table 2 provides a comparison of the computational effectiveness of the two 
approaches for the control of the stiffened laminated composite plate problem for 
cases involving different number of design variables and different number of fre- 
quencies being controlled. 

Table 2. Finite Difference Versus Nelson' Approach - Normalized CPU Time 


Type of 
Design 

Finite Difference 
Approach 

Nelson's Approach 

13 design 
variable, 

3 frequencies 

1.598 

1.0 

13 design 
variables, 

8 frequencies 

1.939 

1.0 

5 design 
variables, 

8 frequencies 

1.700 

1.0 


It is clear from Table 2 that in all the cases considered the finite 
difference approach requires more CPU time as compared to the Nelson's Analytic 
approach. The percentage increase in CPU time increases with the number of 
frequencies considered. In the finite difference approach an eigenvalue problem 
needs to be solved for each design variable considered. If shifting the K matrix 
in conjunction with inverse iteration is used to calculate the 

eigenvalues/eigenvectors, the finite difference approach is quite competitive with 
the Nelson's analytic approach. Even though the approach may require about twice 
the time of Nelson's method, the coding effort is far.,less M in the case of the fi- 
nite difference approach. Secondly, calculation of (derivatives of stiff- 

dp dp 

ness and mass matrices respectively) required in the case of Nelson's approach 
can be quite difficult in some cases... In the case where the design variables are 
element frame areas, calculation fairly straightforward. However, 

if p, corresponds to the number of plies with a given orientation then the calcu- 
lation of is fairly involved. 

In conclusion, it needs to be emphasized that the finite difference scheme 
for the calculation of the eigenvalue and eigenvector derivatives does not appear 
to be costly enough to warrant the use of the analytical method. With the former 
scheme one does not have to "tinker" with the "black box" that generates the 
eigenvalues and eigenvectors for a given design variable vector. The analytical 
method on the other hand needs an intimate knowledge of this "black box". 
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